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Abstract 

We present leading logarithmic QCD corrections to the decay Bg — > 77 in the Standard Model. 
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state effects, and the renormalization scale fj,. 
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1 Introduction 



Rare B decays induced by flavor changing neutral currents (FCNC) are known to provide information 
about the Standard Model (SM) at quantum level and quantitative information on the SM parameters, 
such as the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The CLEO observation of the 
radiative decay mode B Xgj has been analysed in the SM and the rate agrees with the SM-based 
theoretical calculations Another example is Bg — > 07, which is CKM allowed due to the dominant 
CKM matrix element dependence of the decay rate. The calculational procedure of such decay rates 
is to use an effective Hamiltonian obtained by integrating out the top quark and the bosons |^] 

Heff = -AvM E C^{^^)Oi{,l) . (1) 

Here Oi are suitable operators and Cj are Wilson coefficients renormalized at the scale /x. The 
coefficients can be calculated perturbatively. Hadronic matrix elements < y|Oj|S > can be calculated 
using some non-perturbative methods like QCD sum rules, which is one of the powerful methods to 
calculate matrix elements in a model independent way. 

Among rare decays, Bg ^ 77 is a potential candidate to test the SM and search for new physics. 
The final state contains CP-odd and CP-even states, allowing us to study CP violating effects. Mea- 
surement of these odd and even states is a powerful test of the underlying theory, in particular SM. In 
the literature, Bs 77 decay has been investigated earlier in the lowest order and the branching 
ratio is found to be 4.5 • 10~^ in the SM context for = 0.5 GeV and other parameters given in 
Table |, using the constituent quark model. 

In the present work we give the leading logarithmic QCD-improved rates for Bs 77. This 
can be achieved through a matching of the full theory with the effective theory at a scale ^ = niw, 
using the effective Hamiltonian in eq. (|^) , and performing an evolution of the Wilson coefficients from 
m]y down to /i ~ 0(mf,), thus resumming all large logarithms of the form a"^ {mh)log'^ (j^) , where 
m < n (n = 0, 1, 2, . . .). In the leading logarithmic approximation, which we use here, m = n. The 
effective Hamiltonian in eq. (ffl) is identical for 5 — > 57 and for h 577 to this order of Since 
there exists after applying the equations of motion no gauge-invariant FCNC-2-photon operator with 
field dimension < 6, the set of operators (eq. @) is a basis for both decays Q. The bound state 
effects of the Bg meson are modeled through an Heavy Quark Effective Theory (HQET) inspired 
approach following . We estimate further the additional contribution in the decay Bg 77 through 
Bg (pj followed by i?i> — > 7 using Vector Meson Dominance (VMD) Q. In the language of the 
operator basis in eq. @), this contribution involves the operator O7, (see eq. @ below). The decay 

~^ 07 was studied in the literature in the framework of Light-cone QCD sum rules |lO| . We have 
repeated the calculation using the ordinary QCD sum rules including the contribution from the gluon 



condensate. The CP-odd and CP-even amplitudes in Bs 77 are then estimated by considering 
the 0^7 process using a (/)-photon conversion factor supphed by the VMD model. In this part an 
extrapolation from p'"^ = (needed for Bs ^7) to p'"^ = (required for B^ — > 77) is necessary. 
We assume, that the form factor is dominated by a single pole, which is a good approximation for 
light mesons. The decay rate for Bs 77 depends sensitively on the model parameters (m;,, A^) and 
^. For typical values {mb,As) = (5 GeV, 370 MeV) and ^=5 GeV, we get (including long-distance 
effects through O7) the branching ratio 13{Bs ")i))sd+lDoj — l-l^ " 10~^, which is a factor 1.9 
larger compared to the lowest order estimate for the same values of the parameters. However, varying 
(mf,,As) and /x in the allowed range results in significant variation on the branching ratio, yielding 
0.38 • 10-6 < ^(^^ ^ 77)5D+LDo, < 1-43 • 10-^. 

The paper is organized as follows: In section 2 we display the amplitude for Bs ^ 77 in a HQET 
inspired model and present the leading logarithmic QCD corrections. In section 3, we calculate the 
form factor Fi in the decay Bs 07 using QCD sum rules and compare our result with the previous 
result obtained in |]l^] . Section 4 is devoted to the estimate of the Bs ^ 4"^ ^ 77 amplitude in the 
framework of VMD. We discuss the resulting branching ratio B{Bs — > ^^)sd+lDo^ its parametric 
dependence on the model parameters (m^, Ag) and the scale in section 5. 

2 Leading logarithmic improved short-distance contributions in Bg 
77 decay 



The amplitude for the decay Bs — > 77 can be decomposed as |^-|6[ 

A{Bs ^ 77) = e^(A;i)e^(A:2)(A+5M. + iA-ef^.af^kfk^) , (2) 

where the ki and e'({hj) denote the four-momenta and the polarization vectors of the outgoing photons, 
respectively]^. Using the effective Hamiltonian in eq. (|l]), the CP-even {A^) and CP-odd {A^) parts 
in the SM can be written as (for diagrams see fig. |^ and fig. |2|) in a HQET inspired approach |^ : 



'''' '\SAs{mB,-As){ml^f + m's") ' 
77152 



9 J^*^// _|_ : 



eff 



{-mbJ{mb) +msJ{ms))D{n) , 



V27r V3mB^As(mB^ - As) 

- ^Qp{'mq)Cq{lJ^) + -^[mbl\{mb) +msl^{ms))D{fi) \ , (3) 

q 9(m^-'-' -I- nis-'-' ) I 



^We adopt the convention Tr(7f 7'^7"7''75) = 4ief""°'3, with e^^^^ = +1. 

■^In an earher version of this paper the contributions of the operators Oi,3...6 in the irreducible part in and A~ 
were not completely taken into account. This is corrected here. As a second improvement we give the amplitudes in a 
formalism inspired by HQET to estimate the uncertainties coming from the bound state. 



where we have used the unitarity of the CKM-matrix J2i=u ct^is^ib = ^^^d have neglected the 
contribution due to V*gVub *C V^lVtb = A^. In eq. (^) Nc is the colour factor {Nc = 3 for QCD) and 
= I for g = n, c and Qg = — ^ for q = d,s,b. The QCD-corrected Wilson coefficients in leading 
logarithmic approximation |^], Ci,,,^{fi) and Cj^^{fi), enter the amplitudes in the combinations 

Csifi) = Cbif^) = (C3(m) + C4(/x))(A^c + 1) - N.C^ifi) - Ceif^) , 
D{fi) = C^ifi) + Ce{fi)N, . (4) 

While Ci...6(/i) are the coefficients of the operators O1...6, CY^^{fi) is the "effective" coefficient of O7 
and contains renormalization scheme dependent contributions from the four-quark operators O1...6 
in Tieff to the effective vertex in 6 ^ 57. In the NDR scheme, which we use here, Cj^-^{fj,) = 
C-j{ijl) — iC5(/i) — Ce{fi), see |^] for details. The initial values of Ci,,s{fJ-) and Cj^^ (fi) in the SM are 

Ci,3...6{^w) = , 
C2{mw) = 1 , 
C7 M = 4(^31)4 1^^+ 24(x-l)3 ' 

and X = m^/m'^,. For comparison, Ci(mfe) = —0.246, C2{mii) = 1.106, C3{mb) = 0.011, C4(mf,) = 
—0.025, C^{mb) = 0.007, C%{mh) = —0.031 and Cj^^ {niij) = —0.313 for the input values given in 



of Heff 


is given as 


Oi = 


{sLalfJ}La){cLf3l^CLf3), 


02 = 


{sLal^hLI3){cL(jYcLa), 


03 = 


q=u,d,s,c,b 


04 = 


q=u,d,s,c,b 


05 = 


q=u,d,s,c,b 


Oe = 


q=u,d,s^c,b 


O7 = 


167^2 SaCTfiuirribR + 'msL)baF'^'' , 


Os = 


-|^g^2 SaTaf3(yfiv{mbR + msL)hpG' 



afiu 



(6) 



where L and R denote chiral projections, L{R) = 1/2(1 =p 75) and a and (3 are SU{3) colour indices. 
Note that Og does not contribute here in this order of a^. The functions I {nig), J {nig) and A{nig) 
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come from the irreducible diagrams with an internal q type quark propagating, see fig. [l|, and are 
defined as 



2 

TTl 

I{mq) = 1 + ^A(mg) , 
m% — Ami 





Am^ 




Am^ 








2 

A(m„) = - I 2arctan(-^^ ^ ^) - vr I for < 1. (7) 

' ms, J Am^ 

The parameter enters eq. (|3|) through the bound state kinematics. For definiteness, we consider 
the decay Bs = (bs) — > 77. We write the momentum of the 6-quark inside the meson as p = mhv + k, 
where /c is a small residual momentum, v is the 4-velocity, which connects the quark with the meson 
kinematics through P = ms^v and P is the momentum of the meson. In the Bg rest frame, v = 
(1,0,0,0). For the reducible diagrams, see fig. ^, we need to evaluate p.ki and p'.ki, i = 1,2, where 
ki,p' are the momenta of the outgoing photon and s-quark, respectively. Now following |^, we average 
the residual momentum of the 6-quark through 

2mb 

< kakf3> = Y(9al3 " VaVjs) , (8) 

where Ai, A2 are matrix elements from the heavy quark expansion. Using P = p — p' , P.ki = 
y,ki = ^ and the HQET relation g 

ms, =mfe + A, --^(Ai + SAa) (9) 



one gets: 



P-h = ^^^{mB,-K 
p -i^t — 2 « ' 
{ml^^f =p'^ = ml- 3A2 , 



{m^fff^p'^ = {mlfff-ml^+2mBAs- (10) 

The non-perturbative parameter A^ can be related to A, which has been extracted (together with 
Ai) from data on semileptonic B^,B^ decays by and the measured mass difference Am = 



— miB = 90 MeV 1 12], defining A^ = A + Am. The matrix element A2 is well determined from 
the 5^*^^ — mass splitting, A2 = 0.12 GeV^. With the help of eq. the correlated values of A 



and Ai can be transcribed into a correlation between A^g^ and mh. We select 3 representative values 
I {mb,As) = (5.03, 370), (4.91, 480), (4.79, 590) in (GeV,MeV) to study the hadronic uncertainties of 
our approach. Furthermore, we have used the definition 

< 0\s^^^5b\Bs{P) > = ifB^P^., (11) 

which leads together with the off-shellness of the quarks inside the meson to the matrix element of 
the pseudoscalar current 

2 

<Q\-s-i^h\Bs{P)> = -ifB. .r""^- . (12) 

+ ms 

The auxiliary function = g^{rri^^^ ^Kg) is defined as 

= msAml^^ + mlfff + Ks{ml^ - (m^^^ +mf /)2) . (13) 

Note that in the limit rUg, m'^-^J — > rrib^s and using ras^ = mb + mg we recover the result obtained 

by the constituent quark model 0-^, ignoring QCD corrections. Using the above expressions, the 
partial decay width is then given by : 

Now, there are 2 new observations to be made: 
First, the Wilson coefficients in eq. (^) depend on the scale /i. Therefore, since the behaviour of these 
short-distance (SD) coefficients under renormalization is known from the studies of i? — > Xsj 
one can give an improved width for 77 by including the leading logarithmic QCD corrections, 

by renormalizing the coefficients C1...6 and C^^^ from fj, = mw down to the relevant scale fi ~ 0{mb). 
The explicit O^Og) improvement in the decay width T{Bs 77) requires the calculation of a large 
number of virtual corrections, which we have not taken into account. Varying the scale ^ in the 
range ^ < /U < 2mb, one introduces an uncertainty, which can be reduced only when the complete 
next-to-leading order (NLO)-analysis is available, similar to the recently completed calculation for the 
B^Xg^ decay |]. 

The second point concerns the strong dependence of the decay width T{Bg 77) on A^, F ~ ^^(tt) 

s 

in eq. (0) . It originates in the s-quark propagator in the diagram with an intermediate s-quark in fig. ||. 

In the earlier work the authors of e.g. Q evaluated the decay width with mg ~ rriK, assuming that the 

constituent quarks are to be treated as static quarks in the meson. This is a questionable assumption. 

In the HQET inspired approach, this gets replaced by A^, which is well-defined experimentally. This 

formalism implies, that the decay width F(i?rf 77) will involve the parameter A, which avoids the 

unwanted uncertainty on m^. 

''We choose (Ai,A) = (-0.09, 280), (-0.19, 390), (-0.29,500) in (GeV^MeV) from fig. 1 in 0. 

5 



b s 



Figure 1: The generic diagram contributing to 6 — 577 in the effective theory due to the (Fierz 
ordered) four-quark operators. The diagram with interchanged photons is not shown. 




Figure 2: The reducible diagrams contributing to 6 — > S77. The blob denotes the FCNC operator 
O7. The diagrams with interchanged photons are not shown. 

Lowering the scale jj, from = mw to ~ 0{mi,) and enhances the branching ratio B{Bs 77). 
The dependence of the branching ratio as a function of the scale ^ for different values of (m^, A^) is 
discussed in the last section including the Oy-type long-distance (LD) estimate. 

3 QCD sum rule for the Bs — ^ 07 form factor 

3.1 Calculation of the sum rule 

The amplitude for the Bg — ^7 transition A{Bs — (?!>7) =< (t)^\'Heff\Bs > reduces to 

A{B, ^ 07) = e^'Cmb < ^{p')\sa^,Rq%\Bs{p) > (15) 

with the constant C 

^=%ii^*sVt,C^,^^l^). (16) 

where we just take the contribution due to the electromagnetic penguin operator O7 into account and 
put rus = 0, justified by <C mi,. Here e and q are the photon polarization and the (outgoing) 
photon momentum, respectively. Lorentz decomposition gives further: 

+ {etp.q-p^q.et')G{q^) , (17) 
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Parameter 


Value 




1.4 (GeV) 




4.8 (GeV) 


"em 


129 


At 


0.04 




4.09 • 10~i3 (QgV) 




0.2 (GeV) 




5.369 (GeV) 


mt 


175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


. (5) 


0.214 (GeV) 




0.117 


A2 


0.12 (GeV^) 



Table 1: Values of the input parameters used in the numerical calculations unless otherwise specified. 



where p, p' denote the four-momenta of the initial i?s-meson and the outgoing <p, respectively and is 
the polarization vector of the c/i-meson. At the point q'^ = 0, it is enough to calculate -Fi(O), since both 



form factors coincide [13|. Note, that the form factors introduced above are in general functions of 
two variables and p''^. Since cp is on-shell, we abbreviate here and in the following unless otherwise 



stated Fi(g2) == Fi{q^,p 



rrij 



The starting point for the sum rule is the three-point function 



ail 



d'xe'P^-'P'y < 0\T[Ux)T^{0)My)]\0 > , 



(18) 



where Jq = s^aS, J5 = sij^b and = Si^fyfivq^h correspond to the electromagnetic, pseudoscalar 
currents and the penguin operator, respectively. Performing now an operator product expansion 
(OPE) of Ta^, we obtain a perturbative term, the so-called bare loop, and non-perturbative power 
corrections, diagrammatically shown in fig. ^. The bare loop diagram can be obtained using a double 
dispersion relation in p^ and p'^ , 



1 



p{s,s') 



fOO POO 

^""'■^ = ^ L Jo '^''{s - p^){s' - p'2) 



-|- subtractions . 



(19) 



Technically, the spectral density p{s,s') can be calculated by using the Cutkowsky rule, namely, by 
replacing the usual propagator denominator by a delta function: 
■ — > —2iTi6{k'^ — m?)0{kQ). As a result we get 



p{s,s) 



N, 



c 4 



(20) 
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Figure 3: Contributions of perturbation theory and of vacuum condensates to the Bs 07 decay. 
The dashed Hues denote soft gluons. 

OPE enables us further to parametrize the non-perturbative effects in terms of vacuum expectation 
values of gauge-invariant operators up to a certain dimension, the so-called condensates. We consider 
up to dimension-5 operators; i.e. the quark condensate, gluon condensate and the quark-gluon (mixed) 
condensate contributions (fig. ^). This calculation is carried out in the fixed point gauge, i.e. j4^.x^ = 



0. We get 



where 



rr. —n^b - 1 

Tdim-3 = < SS> 



2 (p2 — m'^)p''^ 



a 



el fl — X roo 



Tdim-A = <G > I dx I dy I daa 



JO JO 
_ rrib _ „ . ml ml 

+ 2(p2-m2)2y2] ' (21) 



4 4 

ci = mf,x , 

C2 = mlx'^{l -x-y) , 

C3 = mlx^{3 + y){l-x-y) 

di = mix , 

d2 = x{l - X - y) , 
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ds = y{l-x-y) . (22) 
Here we used the exponential representation for the gluon condensate contribution: 



1 1 



da a"-ie-"^ . (23) 



(n-l)!yo 

The momenta P, P' in eq. (^T|) are euchdean. 

For the calculation of the physical part of the sum rules we insert a complete set of on-shell states 
with the same quantum numbers as and 4> in eq. ([l8|) and get a double dispersion relation 



T,,,. = ^^/0"^0 (^2_^|j(p.2_^2) ^i (0) + continuum , (24) 

where and are the leptonic decay constants of the (j) arid Bg mesons respectively, defined as 
usual by 

< 0| J„|</> > = m^f^ei , 
< Q\MB,{p) > = fBsmljmb . (25) 

We have absorbed all higher order states and resonances in the continuum. 

Now, we equate the hadron-world with the quark-world by Tp^yg = Thare+T^+T^+T^. Using quark- 
hadron duality, we model the continuum contribution by purely perturbative QCD. To be definite, it 
is the part in eq. ( [T9| ) above the so-called continuum thresholds sq and Sq. To get rid of subtractions 
and to suppress the contribution of higher order states, we apply a Double Borel transformation B 
[ p!7| ] with respect to p'^ and p'^. We make use of the following properties of the Borel transform: 

1 (—1)"" e"™^/-^^ 

^((p2_^2)n) = {M^)n ' ^^^^^ 

Bie-'^P') = 6{l-aM^) . (27) 

Finally, this yields the sum rule: 

F.(o) = r- , {4,rd,fd/M«./)e-/M'-vM" 

M'^ jBsUm^m^^^ vr^ J ml Jo 

^ <ss>e [1 "^0(4^4 + 3M2M'2 2M2^J 

+ —<G^> Nix)dx} , (28) 

vr Jo 

where s = min{s — m^, Sg) and x^ax = m'^+M'^ • Here we used the parametrization 

g < saGs > = ml < ss > . (29) 



The last term in eq. ( p^ is due to the gluon condensate contribution and the function N{x) is defined 
by: 

N{x) = 4^ ^^P(- M2(l - x'^ImVM-^) )^lM''-(^lM'^ - ^M'M'^ + 5M'M''x 

+ 5M^M'\-M2M'^x2-2M^M'\2-M^x^)/(M^(-M'2 + M'\ + Af2j;)^). (30) 
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Figure 4: The dependence of the decay constant Fi{0) on the Borel parameters and M'^ for 
so = 33 GeV^. 

3.2 Analysis of the sum rule 

First we hst the values of the input parameters entering the sum rules (eq. (|2^)), which are not 
included in Table |l|: mg = 0.8 GeV^ [|l5|, < ss >= -0.011 GeV^ ll|, < >= 0.03 GeV^ Q, 
= 1.019 GeV and = 0.23 GeV p|. 

We do the calculations for two different continuum threshold values sq = 33 GeV^ and sq = 
35 GeV^ and take Sg = 1.8 GeV^. In fig. |^ we present the dependence of Fi(0) on Af^ and M'^ for 
So = 33 GeV^. According to the QCD sum rules method, it is necessary to find a range of and 
M'^, where the dependence of -Fi(O) on these parameters is very weak and, at the same time, the 
power corrections and the continuum contribution remain under control. From fig. Q and fig. ^ follows 
that the best stability region for Fi(0) is 7 GeV^ < < 9 GeV^ 2 GeV^ < M'^ < 3 GeV^ for 
So = 33,35 GeV^. We get: 

Fi(0) = 0.24 ±0.02 . (31) 
This agrees for our value of rrih within errors with the result given in the literature, based on Light-cone 



QCD sum rule calculations |1C]. 

Numerical analysis shows, as also mentioned in [^], that the natural hierarchy of the bare loop, the 
power corrections and continuum contributions does not hold due to the smallness of the integration 
region, and the power corrections exceed the bare loop contribution. The gluon condensate contribu- 
tion is < 1% of the dim-3 + dim-5 condensate contributions and can therefore be safely neglected in 
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Figure 5: The dependence of the decay constant i^i(O) on the Borel parameter for fixed M'^ at 
So = 33 GeV^ (sohd) and sq = 35 GeV^ (dashed). 

numerical calculations. 

4 The Bg ^ ^ 77 amplitude using VMD model 



We consider the construction of a VMD amplitude using the amplitude for the decay Bg — > 07 as an 
input. Our aim is to continue the Bs — > (p^ decay amplitude from p'^ = to p'"^ = 0, such that 
the (j) meson propagates as a massless virtual particle before converting into a photon. Note that we 
suppressed in our notation the dependence of the form factor Fi{q^) = Fi{q'^,p''^ = m^) on the second 
argument p'"^. We define here Fi{Q'^) = Fi[q^ = 0,Q^) for virtual momenta = —p'"^. Assuming 
pole-type behaviour of the form factor Fi{Q'^) we extrapolate using the single-pole form 

.2^ _ A(0) 



FiiQ' 



l-Q2/m2 



(32) 



pole 



which works well for light mesons. Using an rupoie of order 1.7 — 1.9 GeV, which corresponds to the 
mass of the higher resonances of 0, we estimate Fi{{)) = 0.16 it 0.02. 

With the help of VMD p^] and factorization we can now present the amplitude for Bg 77. 



Using the intermediate propagator 



-1 



at Q2 



0, the 



7 conversion vertex from the VMD 



mechanism 

< 0|J^em|<A(p',e) >= eQsf4,{0)m^efj, , 
and the A{Bs ipj) amplitude, see eq. (p^), we get: 

A{Bs ^ 07 ^ 77) = e^(A;i)e^(fc2)(A+^^^ff^, + i^Z^o, ^2 ) 

with the CP-even {Aj^^^ ) and CP-odd {Aj^j-,^ ) parts: 

2 2 
mf, — m, _ 

Atoo, = 2xCm,^^^F,{0) 
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(33) 



(34) 



= 2V2^^F,{0)Umt^C^7"i^^) . (35) 
vr 3m^ 



where /(/,(0) = 0.18 GeV Q, Qs = —1/3 and C is defined in eq. (16). The factor 2 comes from the 
addition of the diagrams with interchanged photons. Note, that while for the analysis of the sum 
rule for Bg — > (/'7 we have used = /^{m'^), here we take into account the suppression in f(f){Q'^) 
going from = m| to = 0. We treated the polarization vector e"^ as transversal and replaced 
e — > ei, e'f' ^ €2, q ^ ki, p' ^ k2- The conversion factor x is defined as 

X = -eQ.^ . (36) 



Adding this to the short-distance amplitudes (eq. (|3|)), we obtain the Bg —>■ 77 width including 
the 07-type long-distance effects: 



2.2 








, , 1 , 1 1 

(5.03,370) - 

(4.91,480) 


2 








(4.79,590) ' 


1.8 




\ 

\ 

\ 

\ 






1.6 




\ 

\ 

\ 

\ 

\^ 






1.4 




, 1 










4 


6 


8 10 



M (GeV) 

Figure 6: Scale dependence of the ratio i?(/i) defined in eq. (|38|). The solid, short-dashed and long- 
dashed lines correspond to the values {mh,As) in (GeV, MeV) as indicated in the figure. The dotted 
line depicts the suggested choice of the scale from B — > ^^7 studies in NLO ||2|,^. The parameters 
used are given in Table ||. 
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5 Numerical estimates 



(38) 



First we study the leading logarithmic /i-dependence of the ratio 

. _ ^{Bs 7l)ifJ-)sD+LDo^ 

T{Bs ^ -f-f){mw)sD+LDo^ 
In the numerical analysis we neglect the masses of the light quarks. From fig. ^ we find an enhancement 
factor of 1.3 — 2.3 relative to the lowest order result obtained by setting fi = mw, depending on the 
model parameter (m^jA^). Varying fi in the range 2.5 GeV < fi < 10.0 GeV, gives an uncertainty 
AR/R{fi = 5 GeV) ^ ±(17, 19,22)% for A, = (590,480,370) MeV, respectively. Here one can argue, 
that the choice /i = ^ takes into account effectively the bulk of the NLO correction as suggested by 
the NLO calculation for B Xsj [||. 

Table ^ shows the combined fj, and model parameter dependence of the branching ratio 

T{Bs ll)sD+LDo, 



B{Bs ll)sD+LDo, 



(39) 



The dependence of the form factor Fi(m^) on the b quark mass has been extrapolated from fig. 3 [p!o|]. 
Here Fi(0) = 0.14, 0.15, 0.16 has been used for mf, = (5.03, 4.91, 4.79) GeV, respectively. Qualitatively, 



(GeV) 


As = 370 MeV 
nib = 5.03 GeV 


A., = 480 MeV 
nib = 4.91 GeV 


Ks = 590 MeV 
nib = 4.79 GeV 


2.5 
5.0 
10.0 


1.43 • 10^6 
1.18 • 10^6 
0.99 • 10-6 


8.1 • lO"'^ 

6.8- lO"'^ 

5.9- lO"'^ 


5.0 • lO"'^ 
4.3 • lO"'^ 
3.8 • lO"'^ 



Table 2: Branching ratio B{Bs ^^)sd+lDoj fo^' selected values (mb,As) and the renormalization 
scale //. 

the influence of the LD contribution through i?^ — > (/>7 — > 77 reduces the width because of the 
destructive interference of the LD + SD contributions. To quantify this, we define 



B{Bs ^ 77)5D 

with T{Bs — > Jj)sD given in eq. (0). We find, that k lies in the range: 

- 15% <K< -27% , 



(40) 



(41) 



depending mainly on (m;,, A^). 

In conclusion, we have reanalysed the decay rate Bs 'j'j in the SM. We included the lead- 
ing logarithmic QCD corrections and investigated the infiuence of the LD-contributions due to the 
chain i3<j — > 07 ^ 77. Depending on A^, the LD-contributions become sizeable. Other possible LD 
contributions may also arise from the 02-type transitions. 
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The decay rate of Bg 77 depends sensitively on (m^, A^) and fi. Fixing fitofi = ^as suggested 
by the NLO calculation of -B — > Xg^ and varying the model parameter {rrij,, Ag) (see Table we find 
that the branching ratio B{Bs ^^)sd+lDoj is uncertain by a large factor 

0.5 • 10-6 < B{Bs ^ ii)sD+LDo, < 1-4 • 10-6 . (42) 

Improving this requires NLO calculation in the decay rate Bg 77. With the choice of (m^, A^) = 
(5.03 GeV, 370 MeV), the resulting branching ratio (1.4 • 10"^) is substantially larger than what has 



been stated in the literature. The present best limit on the decay Bg 77 is [20| 



B{Bg 77) < 1.48 • 10"^ , (43) 
which is still a factor ~ 100 — 300 away from the estimates given here. 

Note added: Recently, the leading logarithmic QCD corrections for the short-distance part of the 
decay Bg — > 77 have also been calculated by Chang et al. |^l|. They derived the decay rate with the 
full set of operators O1...8 and we agree with their analytical expression. Our model to incorporate the 
bound state effects in the Bg meson is inspired by HQET, resulting in the parameters {mi,, Ag). The 
strong parametric dependence of the decay rate T{Bg 77) on (rn-j,,As) and on fi has been studied 



by us; Chang et al. |21] fix A = m^^ — ruf,, using the naive constituent quark model, and set fi = mi,. 
We emphasize here that the decay rate is sensitive to both of these parameters and requires further 
theoretical investigation. 
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